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Abstract 

The first empirical model of the equatorial mass density of the plasmasphere is 
constructed using ground-based ULF wave diagnostics. Plasmaspheric mass density 
between L=\ .7 and 7=3.2 has been determined using over 5200 hours of data from pairs 
of stations in the MEASURE array of ground magnetometers. The least-squares fit to the 
data as a function of L shows that mass density falls logarithmically with L. Average ion 
mass as a function of L is also estimated by combining the mass density model with 
plasmaspheric electron density profiles determined from the IMAGE Radio Plasma 
Imager (RPI). Additionally, we use the RPI electron density database to examine how 
the average ion mass changes under different levels of geomagnetic activity. We find 
that average ion mass is greatest under the most disturbed conditions. This result 
indicates that heavy ion concentrations are enhanced during large geomagnetic 
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disturbances, and therefore play an important role in storm-time plasmaspheric dynamics. 
The average ion mass is also used to constrain the concentrations of He + and 0 + . 
Estimates of the He + concentration determined this way can be useful for interpreting 
IMAGE Extreme Ultraviolet Imager (EUV) images. 

Introduction 

Many studies have been performed using spacecraft data and VLF whistlers to 
determine the equatorial electron density profile of the inner magnetosphere [Chappell, 

1 972, Farrugia, 1 989, Carpenter and Anderson, 1 992, Fung et al. , 200 1 , Sheeley et al., 
2001]. The usual result is that plasmaspheric equatorial electron density (n<.)eq ranges 
from hundreds to many thousands per cubic centimeter and decreases as a power law 
with radial distance from the Earth. Diurnal [Park, 1978], annual [Clilverd et al, 1991], 
and storm-time [Park, 1973; Park, 1974] variations of the total electron plasma density 
have all been observed. 

The mass density and composition of the plasmasphere are not as well known, 
particularly the inner plasmasphere (L <3). This is because satellites pass through this 
region quickly and spacecraft charging effects make direct measurements of the dense, 
low-energy plasma difficult [Moldwin, 1997]. Most of the mass density and composition 
observations of the plasmasphere have come from the DE 1, GEOS 1 and 2, and ISEE 1 
missions [Horwitz et al., 1984, Farrugia et al., 1989]. These studies have shown that 
He + and 0 + are the major heavy ion species in the plasmasphere and usually have 
concentrations of about 2-6% for He + and 1% for 0 + . Craven et al. [1997] reported 
typical He+ to H+ ratios of ~0.03 to 0.3. Horwitz et al. [1984] found that in the aftermath 
of a storm 0+ density could become comparable to H+ density in the plasmasphere. 
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Techniques for remotely sensing the mass density along closed magnetic field 
lines involve using pairs of ground-based magnetometers to measure field line resonance 
frequencies. These frequencies can then be used to infer the mass density along the 
magnetic field, much like the mass density along a string fixed at both ends can be 
determined from its harmonic frequencies [Schultz, 1996; Denton and Gallagher, 2000 J. 
Two of these techniques are the power ratio technique [Baransky et al., 1985, 1989] and 
the cross-phase technique [Waters etal., 1991, 1994]. Berube etal. [2003] have 
developed an automated method for determining field line eigenfrequencies using pairs 
of ground magnetometers. This new routine allows the analysis of large amounts of data, 
and hence the development of a mass density database for the plasmasphere. 

In this study, a model of the equatorial mass density of the plasmasphere using 
ULF field line eigenfrequencies is constructed using data from pairs of ground-based 
magnetometers. The mass density model is then combined with the IMAGE RPI electron 
density model to determine the average ion mass and estimate the maximum 
concentration of heavy ions in the plasmasphere. 

Equatorial Mass Density Model 

The MEASURE array of ground magnetometers [Moldwin et al, 2004, 
htt p://measure.igpp.ucla.edu] is located along the east coast of the United States and 
consists of ten stations covering L= 1 .7 to L— 3.2. Data from seven of the stations were 
used in this study. Field line eigenfrequencies have been determined using the method of 
Berube et al. [2003] for approximately 5200 hours of data from six pairs of stations in the 
array from the years 1999-2001 (The station pairs straddle £=1.74, 1.99, 2.30, 2.57, 2.89 
and 3.12). The uncertainty associated with the selected frequencies is typically 1-2 mHz, 
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and is discussed in detail by Berube et al. [2003]. Hourly averaged mass density is 
computed from the eigenfrequencies by assuming a r' 3 dependence of density along the 
field line and using the results of Schulz [1996], who provides a relationship between 
field line resonant frequencies and mass density. Field line resonances are primarily 
present only during the daytime due to the reduced conductivity of the nightside 
ionosphere. Therefore, this model is only valid on the dayside. 

Figure 1 shows the hourly mass density estimates obtained for each of the six 
pairs of magnetometers, plotted as a function of I. The line is the least squares linear fit 
to the log average density at each L. It is given by the expression 

_eq(X)= io-°- 67L + 51 1.7 <= L <= 3.2 (I) 

The equation represents the average mass density of the dayside plasmasphere since it 
includes all levels of activity. For quiet times, 

_ eq (L) = 10-°' 65L + 5 '' 1.7 <= L <= 3.2 (2) 

-9 < Dst < -3 nT 

For the most disturbed times (Dst <-100 nT) 

_ eq = 1 0-°' 74L + 5 5 1 .7 <= L <= 3 .2 (3) 

Dst < -100 nT 

Radio Plasma Imager Electron Density 
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The Radio Plasma Imager (RPI) [Reinisch et al., 2001] on board the IMAGE 
satellite makes in situ measurements of magnetospheric electron density. Using all RPI 
equatorial electron density measurements from May, 2000 through May, 2001, Fung et 
al. [2001] have constructed electron density profiles for various conditions and local 
times. The expressions for electron density were determined by first binning the data by 
L, and then performing a least squares fit to the binned data. For all the data, 

ne q = 10 -°- 66L+4 - 89 2 <=/,<= 3.5 (4) 

For quiet times, 

n €q = io"°‘ 51L + 4,56 1.7 <~L<= 3.2 (5) 

-9 < Dst < -3 nT 

and for disturbed times, 

neq= 10' 104L+5 - 41 2<=L<= 3.5 (6) 

Dst < -100 nT 

Average Ion Mass and Composition of the Plasmasphere 

The average ion mass can be determined from the mass and electron density 
profiles since Ma Vg = _eq/neq. Figure 2 is a plot of the average ion mass determined from 
the mass density and electron density models described above. The black line is the 
average ion mass found using all the data. The red and blue lines are the quiet and 
disturbed time profiles, respectively. The lighter weighted lines represent the uncertainty 
associated with each quantity, and are based on the uncertainties in the parameters of the 
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least-squares fit to the data. The overall and quiet time profiles are close to unity for 2 < 
L< 3. On average, the plasmasphere is made up of mostly H + . For the most disturbed 
times (Dst < -100 nT) the average ion mass is significantly larger, and increasing with 
increasing L. During storms, a considerable number of heavy ions are present in the 
plasmasphere, with more heavy ions at larger L, indicating the presence of a heavy ion 
torus in the plasmasphere. Since the major constituents of the plasmasphere are H + , He + , 
and 0 + , their concentrations can be estimated with the following relationship: 

M avg *m H .[%H+} + m He .[%He + ] + m 0 .[%0 + ] (7) 

If only one heavy ion species is assumed to be present, we can determine the maximum 
concentration of that species. For example, the upper limit on the average ion mass at 
L -, 2 overall is approximately 2.2 amu, which implies a maximum He + concentration of 
40% by number and a maximum of 8% for 0 + . Table 1 shows maximum heavy ion 
concentrations determined this way for the entire data set, quiet times, and disturbed 
times. The maximum percentage is determined from the upper limit on average ion mass 
in each case. It should be stressed that this maximum determination only assumes one 
heavy ion species. For example, >100% He + concentration infers that there must be 
heavier ion species present in order to explain the observed average ion mass. 

Discussion 

The He + and 0 + concentrations in Table 1 are estimates based on the average 
mass and electron densities of the plasmasphere. The findings presented here agree well 
with the results of previous studies. For example, Craven et al. [1997], using DEI/RIMS 
data, found He + /H + ratios of -0.03-0.3 (roughly 3-25 percent if we assume no other heavy 
ions are present). They also found that the relative abundance of He + /H + decreases with 
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distance from the Earth. Our findings are fairly consistent with both these results, 
particularly for the quiet-time plasmasphere. Recently, Dent et al. [2003], using the same 
ground techniques as this study combined with simultaneous electron density 
measurements, for a single case study, found quiet-time He + concentrations of 35-64 
percent and 0 + concentrations of 7-13 percent for L < 3.45, and almost no heavy ions for 
L > 3.45. Our average results are slightly lower than their single event estimates, but the 
trend of less heavy ions at higher L for quiet times is the same. 

The increase in heavy ion concentrations during disturbed times has interesting 
implications for ring current and radiation belt dynamics. In particular, the results show 
that during severe storms ( Dst< -100 nT), the plasmasphere is dominated by heavy ions, 
rather than H + . 

Finally, independent estimates of the abundance of He + in the magnetosphere are 
a necessary aid to the inversion of IMAGE Extreme Ultraviolet (EUV) images. 
Information about the mass density along the line-of-sight is necessary to accurately 
determine the He + density distribution in the plasmasphere. Clilverd et al. [2003] found 
for two events that changes in the He + column abundance appear to correlate with 
changes in He + determined from ground-based VLF and ULF observations. Future 
studies will combine ground-based ULF wave computations of mass density with in situ 
and VLF-derived electron densities and IMAGE EUV He + observations to create a more 
complete picture of the plasmasphere. 
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MEASURE hourly averaged mass density 



Figure 1. Average plasmaspheric mass density from 6 pairs of stations in the 
MEASURE array. The line represents the least squares linear fit to the mean density at 
each L and is given by equation 1 in the text. 



12 


Average Ion Mass 1.7 < L < 3.1 



Figure 2. Average ion mass derived from equatorial mass and electron density models. 
The three lines represent the ion mass overall, low activity (-9 nT < Dst < -3 nT), and 
disturbed activity (Dst , <-100 nT). The lighter dashed lines represent the uncertainty in 
each quantity, based on the uncertainties in the parameters of the least squares fit to the 
data. 
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He + 

0 + 


L=2 

L=3 

L=2 

L=3 

All Data 

40% 

45% 

8.0% 

8.9% 

Quiet Times 

43% 

29% 

8.6% 

5.7% 

Disturbed 

Times 

100% 

100% 

25% 

62% 


Table 1. Maximum heavy ion concentrations determined from ground-data-derived mass 
density model and RPI in situ electron density profiles, using equation (7) from the text, 
and assuming only one heavy ion species is present for each case. 













